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Abstract — In this paper we introduce the notion of optimiza- 
tion under control and communication constraint in a robotic 
network. Starting from a general setup, we focus our attention 
on the problem of achieving rendezvous in minimum time 
for a network of first order agents with bounded inputs and 
limited range communication. We propose two dynamic control 
and communication laws. These laws are based on consensus 
algorithms for distributed computation of the minimal enclosing 
ball and orthotope of a set of points. We prove that these control 
laws converge to the optimal solution of the centralized problem 
(i.e., when no communication constrains are enforced) as the 
bound on the control input goes to zero. Moreover, we give a 
bound for the time complexity of one of the two laws. 

I. Introduction 

The interesting aspect of motion coordination consists in 
combining together problems from control and communica- 
tion theory. The main difficulty deals with integrating the 
sensing, computing, communication and control aspects of 
problems involving groups of mobile agents. A well known 
problem in control theory is optimal control. Roughly speak- 
ing, it consists in finding a feedback law that minimizes some 
cost functional under some inputs and dynamics constraint. 
In this paper we introduce the notion of optimal control and 
communication for a network of robotic agents. We want 
to study how to solve an optimization problem, in presence 
of both the usual motion constraints and the communication 
ones. In particular this paper is a preliminary contribution 
towards what might be loosely referred to as "distributed 
geometric optimization." In fact many optimization problems 
for robotic networks can be shown to be equivalent to the 
computation of geometric shapes. While in a centralized 
setting the solution is usually simple, the problem becomes 
very complicated when it must be solved in a distributed 
way. Distributed computation over network has been largely 
studied for fixed topologies; e.g., see [1]. 

In this paper we point our attention on the well known 
rendezvous coordination task and look for solutions that 
solve such task in minimum time. We look for distributed 
solutions in networks of mobile agents with first order 
dynamics, bounded inputs and limited-range communication. 

The "multi-agent rendezvous" problem and a "circumcen- 
ter algorithm" have been introduced by Ando and coworkers 
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in [2]. The algorithm proposed in [2] has been extended to 
various synchronous and asynchronous stop-and-go strate- 
gies in [3]. A related algorithm, in which connectivity con- 
straints are not imposed, is proposed in [4]. In [5] the class 
of "circumcenter algorithms" has been studied in networks 
of agents whose state space is R d , for arbitrary d, and with 
communication topology characterized by proximity graphs 
spatially distributed over the disk graph. In [6] the (time and 
communication) complexity of this and other algorithms has 
been studied. All these coordination schemes are memoryless 
(static feedback). In this paper we want to explore dynamic 
control and communication laws in order to approximate 
the optimal solution of the minimum time rendezvous. In 
particular the control and communication laws is based 
on reaching consensus on some logic variables and at the 
same time moving toward the current estimation. A similar 
approach was used in [7] where the agents try to reach a 
consensus on a set of variables called coordination variables. 

Studying the minimum time rendezvous problem in the 
centralized setting we show that, depending on the norm 
used to bound the control input, the optimal solution consists 
of moving toward the center of the minimal enclosing ball 
(bound on L2 norm) and toward the center of the minimal 
enclosing orthotope (bound on the infinity norm) of the 
points located at the initial position of the agents. 

Our main result is the design of a control and com- 
munication law based on a consensus algorithm for the 
distributed computation of the minimal enclosing ball and 
the minimal enclosing orthotope of a set of points. We prove 
the correctness of the two consensus algorithms and provide 
a bound on the time of convergence for the orthotope case. 
Then we prove that the control and communication law that 
combines the consensus with the motion law converges to the 
optimal solution as the control bound goes to zero. Moreover, 
for the problem with input bounded by the infinity norm 
(corresponding to the computation of the minimal enclosing 
orthotope), we prove that the control and communication law 
is a constant factor approximation of the centralized optimal 
solution. 

In Section [H] we introduce a formal model of robotic 
network inspired by the one introduced in [6]. Moreover, 
we define the optimal control and communication problem. 
In Section m] we characterize the solution of minimum time 
rendezvous in a centralized setting. In Section HV1 we define 
the FloodMEB and FloodMEO algorithms for the distributed 
computation of minimal enclosing ball and orthotope, prove 
their correctness and give bounds on time complexity. Sec- 
tion|V]contains the control and communication laws based on 
the consensus algorithms described in Section HVl Finally in 



Section fV-CI and Section IVT1 we show simulations and draw 
the conclusions with future perspectives. 

II. Preliminary developments 

In this section we recall the concepts of network of 
robotic agents, coordination tasks and complexity measures, 
and introduce the notion of optimization under motion and 
communication constraints. 

A. Notation 

We let N, No, and R + denote the natural numbers, the 
non-negative integer numbers, and the positive real numbers, 
respectively. We let Hie{i «} &i denote the Cartesian 
product of sets Si, . ■ ■ , S n . For p £ R, we let \p\ and \p] 
denote the floor and ceil of p. For r £ R+ and p £ R d , 
we let B(p, r) denote the closed ball centered at p with 
radius r, i.e., B(p,r) — {q £ M. d | \\p — q\\2 < r} and 
C(p,r) denote the closed hypercube centered at p with 
sides of length r and parallel to the coordinate axes, i.e., 
C(p,r) = {qeR d \\\p-q\\ oc <r}. 

For /, g : N — ► R, we say that / £ 0(g) (respectively, 
/ £ if there exist no £ N and k £ M + such that 

|/(n)| < k\g(n)\ for all n > uq (respectively, \f(n)\ > 
k\g(n)\ for all n > no). If / £ 0(g) and / £ £l(g), then we 
use the notation / £ &(g). 

Next, we briefly review some useful proximity graphs. 
Given r cmm £ R + , the disk graph <7disk(rcmm), respec- 
tively cube graph £ C ube(7cmm), is the state dependent 
graph on W l defined by the following statement: for any 
pointset {pW, . . . ,p'™'} C K d , the pair (i,j) is an edge 

) • ({p [l1 , . . . ,p [n] }), respectively £ C ube(V cmm J 
({pW, . . . ,p^ n ty), if and only if i ^ j and 

|| p W _ p U] || 2 < r cmm ^ pM - pbl £ S(0 d , r cmm ), 

respectively 

||pW _ p \j] 1^ < rcmm p W _ p W e C (0 d , r cmm ). 

Another useful graph is the complete graph Q cmp \, i.e., the 
graph with edges between any pair of nodes. 

Finally, given a graph Q (even not state dependent), we 
denote with distg(j,j) the topological distance between i 
and j, i.e., the minimum number of agents to go from i to j 
in the graph Q. We define diamg, the diameter of Q, to be 
the maximum topological distance, distg (i,j), for all (i,j). 

B. Modeling a network of robotic agents 

We describe a (uniform) network of robotic agents using 
the formal model introduced in [6] modified for the discrete 
time case. The network is modeled as a tuple (I,A,E cmm ). 
I = {l,...,n} is the set of unique identifiers (UIDs); 
A = {A®} ie i = {(X,U,X ,f)} ieI is called the set of 
physical agents and is a set of control systems consisting of 
a differentiable manifold X (state space), a compact subset 
U of R m (input space), a subset X$ of X (set of allowable 
initial states) and a (sufficiently smooth) map f : Xy.U — > X 
describing the dynamics of ith agent; E cmm : X n — > / x I 
is called the communication edge map. 

The robotic network evolves according to a discrete-time 
communication and motion model. 



Definition 2.1 (Control and communication law): Let S 
be a robotic network. A (uniform, synchronous, dynamic) 
control and communication law CC for S consists of the 
sets: 

(i) L, a set containing the null element, called the 
communication language; elements of L are called 
messages; 

(ii) W, set of values of some logic variables w^ 1 ', i £ 7; 

(iii) Wq C W, subsets of allowable initial values; 
and of the maps: 

(i) msg : X x W x / — > L, message-generation function; 

(ii) stf : W x L n — > VF, called state-transition function; 

(iii) ctl : X x W x L" — > U, called control function. □ 
Roughly speaking this definition has the following mean- 
ing: for all i £ I, to the ith physical agent corresponds a 
logic process, labeled i, that performs the following actions. 
First, at each communication round the ith logic process 
sends to each of its neighbors in the communication graph a 
message (possibly the null message) computed by applying 
the message-generation function to the current values of 
and w^. After a negligible period of time, the ith logic 
process resets the value of its logic variables by applying 
the state-transition function to the current value of w^K and 
to the messages received at time t. Between communication 
instants, the motion of the zth agent is determined by 
applying the control function to the current value of and 
the current value of vM. This idea is formalized as follows. 

Definition 2.2 (Evolution of a robotic network): Let S be 
a robotic network and CC be a control and communication 
law for S. The evolution of (5, CC) from initial conditions 
x £ Xq and w £ Wo, I £ I, is the set of curves x m : 
N — > X and w$ : N — > W, i £ I, satisfying 

x®(t + l) = f(x®(t), ctl(x® (t),w®(t+ l), y W (i))), 

where, for i £ 7, 

+ = stf(u>M 

with the conventions that x^(to) = Xq and wM(io) = Wq, 
i £ I. Here, the function j/M : N — > L™ (describing the 
messages received by agent i) has components 

y [i\, t , = {msg(xW(t),wW(t),i), if (i,j) GE cmm , Q 
J [null, otherwise. 

In the paper we consider the following network. Each 
agent i occupies a location p 1 ^ £ R d , d £ N, and moves 
according to the first order discrete-time integrator 

pM(t+l)=pW(t)+uM(t). (1) 

The communication edge map can be either the one arising 
according to the disk graph, Edisk, or the one according 
to the cube graph, E cu b e . Each control takes values in 
a bounded subset of R rf , that can be either B(0, r clr ) or 
C(0,r ctr ), i.e., ||«M|| 2 < r ctr or < r ctr . Notice that, 

in general, the type of communication edge map and the 
type of control bound are not related. Finally the control 
and communication law will be defined depending on the 
coordination task. 



C. Coordination tasks and time complexity 

We are ready to define the notion of task and of task 
achievement by a robotic network. 

Definition 2.3 ( Coordination task): Let S be a robotic 
network. A (static) coordination task for S is a map T : 
X n — > {true, false}. Additionally, let CC a control and 
communication law for S. The law CC achieves the task Tif, 
for all initial conditions x$ G X and w$ G Wq, i £ I, the 
corresponding network evolution t i— > (x(t),w(t)) has the 
property that there exists T G N such that T(x(t)) = true 
for all t > T. □ 

We are finally ready to define the notion of time com- 
plexity as the minimum number of communication rounds 
needed by the agents to achieve the task Twith CC. 

Definition 2.4 (Time complexity): Let S be a robotic net- 
work and let Tbe a coordination task for S. Let CC be a 
control and communication law for S compatible with T. 
The time complexity to achieve T with CC from xo G Xq is 

TC(T, CC,x„) = inf {T e N \ T{x(t)) = true , V t > T} 

where 1 1— ► (x(t),w(t)) is the evolution of (S,CC) from the 
initial condition (xo,wq). 

The time complexity to achieve Twith CC, TC(T, CC), is 
the maximum TC(T, CC,xo) over all initial conditions xq. 

D. Optimal control and communication in robotic networks 

Having defined a coordination task for a robotic network, 
we can ask whether such task can be accomplished minimiz- 
ing some cost functional. In what follows we will introduce 
the notion of optimal control and communication problem 
and of optimal control and communication law as solution 
of the problem. 

Definition 2.5 (Optimal control and communication): 
Given a task T and a cost functional J(u(-),x(T),T), an 
optimal control and communication problem is the following: 

mrrdrnize u( .) i;I .(o) iX (T),T J(u(-),x(T),T) 

J(«(-), x{T),T) = £^ =0 (/(x(t), u (r)) + g(x(T)), 

subj. to 

(i) (x(-), it(-)) is an input-state trajectory of A, 

A={A®} ieI ; 

(ii) i and j can communicate if and only if 

(i,j)€E cmra (xW(t),...,xW(t)); 

(iii) T(x(t)) = true for alii > T, T G N. 

where I : X n x U n — > K is a sufficiently smooth and 
nonnegative-valued function, called stage cost, and g : 
X n — » R has the same properties plus g(x) = for all 
x G X n such that T(x) = true (for an admissible CC). □ 

We say that a control and communication law CC is 
optimal with respect to the coordination task T and the 
cost functional J, if it solves the above optimal control and 
communication problem. 

We call CC a centralized optimal control and communica- 
tion law if it solves the optimization problem for a network of 
robotic agents that communicate according to the complete 
graph, i.e., the communication edge map is E cmp \. 



Remark 2.6: The centralized solution of an optimal con- 
trol and communication problem is the classical solution of 
the optimal control problem for the whole network system 
without communication constraints. □ 

III. Centralized minimum time rendezvous 

In this section we study the rendezvous problem for a 
robotic network of first order agents with communication 
edge map Edisk or -E cu b e and look for a control and commu- 
nication law that solves the problem in minimum time. 

More formally, let S = (J, A, E cmm ) be a uniform robotic 
network. The (exact) rendezvous task T m d ZYS : X n — > 
{true, false} for S is the static task defined by 

{true, if x® = x®, 
V(i,j) G E cmm (x), 
false, otherwise. 

forx= (x [1 \...,x [n] ). 

Thus, given the uniform network S = (I,A,E cmm ), the 
minimum time rendezvous problem for first order agents 
with limited-range communication and bounded control 
input is the following: 

minimize J2t=o ^ 
u(-),p(T) 
subj. to 

(i) (p(-),u(-)) is an input-state trajectory of A, 

A = {AW}< e i = {(R d , U, R d , /)} ieJ , p(0) = p ; 

(ii) i and j can communicate if and only if 
(i,j)£E cmm (j>m(t),..., P W(t)); 

(iii) T adzvs (pW,. . . = true for all t > T, T G N. 
Here U is either 5(0,^) or C(0,r ctr ), f(p [ii (t), ufi(t)) = 

+vfi(i) and the communication edge map E cmm is 
either £d isk or E cube . 

We refer to the minimum time rendezvous problem 
with communication edge map E cmm and input set U as 
MTK(E cmmi U). 

Next, we provide some preliminary results for 
the centralized setting of the above problem, that 
is, for MTTZ{E cm?h U). Let MEB(pW • • -pW) and 
MEO(pW • • - pi"!) the minimal enclosing ball and orthotope 
of points (pW-.-pW), and let MBC(pW • • -pW) and 
MOC(pW---pN) the centers of MEB(pW • ■ -pW) and 
MEO(p[ 1 ' • • -p[ n l) respectively. We present the following 
theorem omitting the proof based on geometric arguments 
because of space constraints. 

Theorem 3.1: For all r ctr G K+, Pq 1 G W l , i G {1, . . . , n} 
the solution of MTK(E cmph U), U = B(Q,r ctr ) or U = 
O(0, r c tr)> is not unique (the problem is not normal). If 
G S(0,r ctr ), i G {l,...,n}, then 

(i) p(T) = p rndzvs , disk - MBC(pW(0), . . . ,p[' l l(0)), 

U M(f) =min{r ctr ,||p rndzvs -pH(i)|| 2 } 

• vers(p rndzvs -p w (t)), i G {1, . . . , n}, 

is a solution of MT72.(£? cmp i,5(0,r c tr)); 



(ii) if Vi G {l,...,n}, \\p® - MBC(pM . . - p M)|| 2 < 
r ctr , then the solution of MT7Z(E cmp i, B(0, r ctr )) 
is given by p(T) = p mAzvs , disk, Pmdzvs, disk 6 

n ie {i „}B(pM,r ctr ), and uM(i) = p rndz vs, disk - 

pM(t). 

Alternatively, if G C(0, r ctr ), ie{l,...,n}, then 

(iii) P{T) = Pmdzvs, cube, Prndzvs, cube G T^MT, M ''' (i) = 

min{r ctr , |p m dzvs,a - Pa W (*) | } sign(p rndz vs,a - p a ^(t)), 
i G {1, ...,n}, a G {1, . . . , d} is a solution of 
WTR(V,C(0,r tt )), where 

MOC(pW(0),...,pW(0))-i(W-i.), 



^-mt =n 

a 

MOC(p [1] (0), 



,P W (0)) + ^(Ux-Za) 



Lax is the largest side of MEO(pW(0), . . . ,p [n] (0)) 
and l d is the side in direction a; 
(iv) if Vi G {l,...,n} ||pW -MOC(pW •••p [ " 1 )||oo < r ctr 
the solution of MTlZ(E cm? i, C(0, r ctr )) is given by 
p(T) = pmdzvs, Pmdzvs G n ie{ i „} C (p M , r ctr ), and 

U M (t) =Pvndzvs-p [il W. □ 

IV. Distributed consensus on minimal enclosing 

BALL AND ORTHOTOPE 

In the previous section we have shown that minimal 
enclosing shapes play a key role in the solution of minimum 
time rendezvous. In fact if the agents could know the center 
of such shapes (ball or orthotope) the solution of minimum 
time rendezvous would be just a control law that drives each 
agent to this point. Therefore, in this section, we want to 
explore two consensus algorithms to compute the minimal 
enclosing ball and the minimal enclosing orthotope of a set 
of given points in K. d in a distributed way. 

Here is an informal description of what we shall refer to 
as the FloodMEB algorithm: 

[Informal description] Each agent initializes the 
minimal enclosing ball to its initial position, then, 
at each communication round, performs the fol- 
lowing tasks: (i) it acquires from its neighbors 
(a message represented by) the coordinates of the 
minimum set of points describing the boundary of 
their minimal enclosing ball and the coordinates of 
their initial position; (ii) it computes the minimal 
enclosing ball of the point set comprised of its and 
its neighbors' set of points and its initial position 
(that it maintains in memory); (iii) it updates its 
logic variables and message as in (i). 
Before describing the algorithm more formally, we need 
to introduce some notation and state some properties of 
the minimal enclosing ball. Given a set of m points 
{qi,...,q m } G R d in generic positions, we denote with 
MEBbndry({gi, • • • > Qm}) the minimum set of points on the 
boundary of MEB({gi, . . . , q m }) that uniquely identify such 
boundary. When the points are in generic position, we let 
MEBbndry({9i, • ■ ■ , Qm}) denote a minimum set of points 
on the boundary of MEB({gi, . . . , Q m }) that identify such 
boundary. Moreover, let MBR : F(R d ) — ► R the function 



that associates to a set of points the radius of the minimal 
enclosing ball of such points. 

Lemma 4.1 (MEB properties): Let Q n a set of n points. 
The following statements hold. 

(i) there exists a subset Qd C Q n of d + 1 elements such 
that MEB(Q d ) = MEB(Q„); 

(ii) for all Q ni ,Qn 2 C Q n with Q ni C Q n2 , then 
MBR(Q ni ) < MBR(Q„ 2 ); 

(iii) if MBR(Q ni ) = MBR(Q„ 2 ), then MBC(Q ni ) = 
MBC(Q„ 2 ); 

(iv) the number of possible values of MBR(Q ni ), for all 
Qni G Q n , is finite. □ 

Remark 4.2: An important implication of Lemma l4TTT i) is 
that MEBbndry({<Zi, ■ ■ ■ , Qn}) has at most d + 1 points, then 
the number of packets in the message sent and stored by 
each agent is at most d + 1 and does not depend on n. □ 

The algorithm is described formally in the following table. 



Name: 


FloodMEB algorithm. 


Goal: 


Solve the problem of computing min- 




imal enclosing ball of a set of points. 


Logic state: 


= (PWbndry,^ 1 ); 


Msg function: msg(x^ , toW , i) = u>M ; 


Initialization: P w b nd ry (0) = {p H (0)}, 




pW(0)=pH(0). 


For i G {1, . . . 


n}, agent i executes at each time t G N: 


1: acquire w^-- 




2: compute 




^bndry(^ 


f 1) = MEB bndry (WV (i) (i)), 


= 


= {wW(t) j G Af(i) U i} 


3: update w^ 1 





Remark 4.3: For the algorithm to converge it is important 
that each agent keeps in memory the coordinates of its initial 
position and thus computes the minimal enclosing ball on the 
points received from its neighbors together with the point 
located in its initial position. In fact a point on the boundary 
on the minimal enclosing ball of n% points is not ensured to 
be on the boundary of the ball of n-i < Hi points. This means 
that the coordinates of the agents on the boundary could be 
taken out from the logic variables during the first iterations. 
This does not happen, for example, for the minimal enclosing 
orthotope. The result is a simplified consensus algorithm. □ 

We are now ready to prove the algorithm's correctness. 

Theorem 4.4 (Correct MEB computation): Let S be a 
robotic network such that the agents can communicate ac- 
cording to some communication edge map E cmm . For any CC 
such that the graph remains connected along the evolution, 
the FloodMEB algorithm achieves consensus on minimal 
enclosing ball. □ 
Proof: In order to prove correctness of the algorithm, 
observe, first of all, that each law at every node converges 
in a finite number of steps. In fact, using Lemma |4~TI each 



sequence corresponds to a ball whose radius is monotone 
nondecreasing, upper bounded and can assume a finite num- 
ber of values. Then we proceed by contradiction to prove 
that all the laws converge to the same ball (same radius and 
center) and that it is exactly the minimal enclosing ball of the 
n points. Suppose that the algorithm converges on different 
balls (different radius or different center) for different agents. 
Then there must exist two agents that are neighbors and have 
different logic variables (corresponding to different balls). 
But this means that, at the following time instant, they have 
to compute the minimal enclosing ball of a larger set of 
points, then either one of them will take the value of the 
other or both of them will change their value and take a 
common one. Iterating this argument we obtain that all the 
agents must converge to a common value. Now, the ball at 
each node contains, by construction, the initial position of 
that node. Since the ball is the same for each node, it contains 
all the initial positions, then it is the minimal enclosing ball 
of the initial positions. ■ 
Remark 4.5: If we admit that the agents have different 
priority, the initial positions of the agents can be shared by 
all the agents in time of order 0(n 2 ). The algorithm is the 
following. Each agent sends the position of the agent with 
higher (or equivalently lower) priority that he has in memory. 
Each position takes time O(n) to spread in the network, 
therefore the total time complexity is 0(n 2 ). Even if we did 
not provide any bound for the time complexity of FloodMEB 
algorithm, however simulations suggest that it should be of 
order Q(n). Moreover, while the algorithm for sharing the 
initial position needs to store a number of packets of order 
O(n), the FloodMEB algorithm needs to store only d + 2 
packets. □ 
Here is an informal description of what we shall refer to 
as the FloodMEO algorithm: 

[Informal description] Each agent initializes the 
minimal enclosing orthotope to its initial position, 
then, at each communication round, performs the 
following tasks: (i) it acquires from its neigh- 
bors a message represented by the coordinates of 
their current minimal enclosing orthotope; (ii) it 
computes the minimal of its and its neighbors' 
enclosing orthotopes; (iii) it stores as new message 
the coordinates of the minimal enclosing orthotope 
computed at the previous step. 
A more formal description of the algorithm is given in the 
following table. 

In the following theorem we prove the correctness of this 
algorithm, together with the fact that it reaches consensus in 
minimum time. 

Theorem 4.6 (Correct MEO computation): Let S be a 
robotic network such that the agents can communicate ac- 
cording to some communication edge map E cmm . For any CC 
such that the graph remains connected along the evolution, 
then the FloodMEO algorithm achieves consensus on mini- 
mal enclosing orthotope. Moreover, it achieves consensus in 
minimum number of communication rounds given by 

T F ioodMEO= max max {dist g (i, i min>a ), dist g (i, i max>a )} 

ae{l,...,d} ie{l,...,n} 



Name 
Goal: 



FloodMEO algorithm. 



Solve the problem of computing min- 
imal enclosing orthotope of a set of 
points. 

Logic state: w w = {wS ] } a &{i,...,d} 

ci [i] [i] \i 

— \\Pmin,a^Pm.a^,ajSa^{l,...,d} 

Msg function: msg(xW , vM , i) = w<- % ' 

Initialization: p^ n>o (0) = (0), 
pL,»(0) =Pa(0) 



For i £ {1, . . . , n}, agent i executes at each time 

1: acquire w"', j E N(i) 
2: compute Va E {1, . . . , d} 

l»L,a( ! + 1 ) = ™%W«U{i}{PmL, a (*)} 
Pmax,a(t+ 1) = maXj- eA r(i)U{i} {Pmax,a(t) } 

3: update Va e {1, . . . , d) 

Wa i Ht + i) = (p [ l n jt + i),p% ax , a (t + i)) 



are the agents that characterize the 
boundary of the orthotope in direction a and minimize the 
topological distance from i. 

The time complexity of the algorithm is of order 0(n).D 
Proof: In order to prove the correctness and the time 
complexity of the algorithm described in Table ??, we need 
to prove that it is equivalent to 2d FloodMax algorithms for 
leader election (two for each direction) running simultane- 
ously. Once we have proven that, the results on correctness 
and time complexity follow from Chapter 4 in [1]. 

The algorithm is clearly a set of FloodMax algorithms for 
leader election. In fact the boundary of the orthotope in each 
direction a is given by the coordinates of the points on such 
boundary which are characterized by the property of having 
the maximum and minimum value of the ath coordinate 
respectively. 

In order to prove that the exact number of communica- 
tion rounds needed is Tfi ooc /MEO, simply observe that it is 
exactly the minimum time for all the leaders to propagate 
their information through all the network. Hence this is the 
minimum time for every possible consensus algorithm to 
converge. But this is exactly the time taken by 2d FloodMax 
algorithms running simultaneously and therefore the time 
taken by FloodMEO. ■ 

In the following lemma we give, for both FloodMEB and 
FloodMEO algorithms, a bound on the time needed by each 
agent to decide that the algorithm has reached consensus. 

Lemma 4.7 (Termination condition): Consider a network 
S, where the FloodMEB (FloodMEO) algorithm is running. 
Each agent can decide that the algorithm has reached con- 
sensus if the value of its MEB (MEO) has not changed after 
diamg communication rounds. □ 
Proof: In order to prove the claim we proceed by con- 
tradiction. Suppose that after diamg communication rounds 



the MEB (MEO) of agent i (for some i E {1, . . . , n}) has 
not changed and the algorithm has not converged yet. Then 
there will exist a T > diamg; such that the MEB (MEO) of 
agent i will change to a new value. But this means that the 
new value, stored T rounds before by some other agent j, 
took a number of communication rounds greater than diamg; 
to arrive from j to i and this contradicts the definition of 
diameter of Q. ■ 

V. Constant factor approximation of minimum 

TIME RENDEZVOUS CONTROL AND COMMUNICATION LAW 

The centralize solution for minimum time rendezvous and 
the consensus algorithms studied in the previous section 
suggest a dynamic control and communication law that plays 
a key role in the minimum time rendezvous problem. 

Here is an informal description of what we shall refer to as 
the move-toward-MSC (MOC) control and communication 
law, CCmeb (CCmeo): 

Each agent initializes its logic variables to its 
initial position, then, at each communication round, 
performs the following tasks: (i) it acquires from its 
neighbors a message given by their logic variables 
and positions; (ii) it runs, as state transition func- 
tion, the FloodMEB(MEO) algorithm; (iii) it moves 
toward the center of the current ball (orthotope) 
while maintaining connectivity. 
Next, we formally define the law as follows. First we 
assume that each agent operates according to the stan- 
dard message-generation function, that is msg(x^ , , i) = 
(x®,w®). Second, before the FloodMEB (or FloodMEO) 
algorithm reach consensus, connectivity is maintained by 
restricting the allowable motion of each agent in some 
appropriate manner. The exact algorithm can be found for 
example in [6]. 

The state transition function implements the FloodMEB 
and FloodMEO algorithms respectively, with logic variables 
as defined in the two tables above. 



Define the control function ctl 
each agent i E {!,... ,n} by: 



x U 



for 



ctl(pM W) = max{A* • (p r ndzvs(w® , y W ) - p W ), r ctr } 

■ vcrs(fw Z f s {w [l] , y W ) - p [l] ), 

with 

Pmdzv S (w [l \y [t] ) = MBC( W W,2/H) (2) 

and A* is chosen in order to maintain connectivity until 
consensus is reached. 

In a network with communication edge map E cmm = E cu \, e 
the procedure described above is applied separately in every 
direction a G {1, . . . , d}. 

The correctness of the two control and communication 
laws is proven in the following lemma. 

Lemma 5.1 (Correctness of CCmeb and CCmeo)- On the 
network S with communication edge map Pbaii or E cu \, e 
and bound on the ith control input G P(0,r ctr ) or 
G C(0, r ctr ), the control and communication laws CCmeb 
and CCmeo achieve rendezvous at MBC(pM (0), . . . , (0)) 
and MOC(pW(0), . . . ,p [n] (0)) respectively. □ 



Proof: By the connectivity arguments done before 
and by Theorem 14.41 and Theorem 14.61 we know that there 
exists T G N such that for t = T the network is 
connected and all the agents have reached consensus on 
MBC(pW(0),...,pN(o)) (or MOC(pW(0), . . . ,P W (())))■ 
Since this instant all the agents can move toward the same 
point (at maximum speed) without enforcing connectivity 
constraint anymore. Thus, they can converge to the ren- 
dezvous point which is exactly MBC(pW(0), . . . ,JJ^(0)) 
(or MOC(pW(0),...,pW(0))). ■ 

A. Time complexity of CCmeb ond CCmeo 

In the previous lemma we have proven that the control and 
communication laws CCmeb and CCmeo achieve consensus. 
Now we ask how fast these laws are depending on the control 
bound r ctr and the number of agents. 

Theorem 5.2: For r cmm E K + , d G N, consider the 
network S with communication edge map either Pdisk or 
-Ecube- The following statements hold: 

(i) for G B(0, r ctr ), i G {1, ...,»}, the control and 
communication law CCmeb asymptotically converges 
to the minimum time rendezvous centralized solution 
MTK(E cmph P(0, r^)) as r ct[ -> 0+ (for all fixed n). 

(ii) for G C(0,r ctr ), i G {1, ...,n}, the con- 
trol and communication law CCmeo converges to 
the minimum time rendezvous centralized solution 
MTTZ(E cmph C(0,r ctI )) for r ctt -> 0+ (for all fixed 
n). Moreover, it is a constant factor approximation 
of MTn(E cmph C(0, r ctr )), i.e., TC(7; ndzvs , CC ME o) G 
9( JL) for r ctI -> 0+ and n -► +oo. □ 

Before proving the theorem let us state a useful lemma. 

Lemma 5.3 ([8]): For all pointsets Pi C P, we have 
MBC(Pi) G MEB(P) and MOC(Pi) G MEO(P). □ 
Proof: [Theorem 15.21 The line of proof of the two 
statements is the same, hence we prove only the second 
one which is stronger due to the stronger result on the time 
complexity of the FloodMEO algorithm. Using the previous 
lemma we know that for all t G No, then P^ndzvs^) e 
MEO(pW(0),...,p[™](0)), where pf ndzvs is defined as in 
This implies that, once the consensus is reached, the 
time to rendezvous is upper bounded by the time of the 
centralized solution. Hence the following bound on the time 
of convergence of CCmeo, Pmeo holds: 



Tj 



MEO 



< 



diam(p [11 (0),...,p [n] (0)) 



TfIoocIMEO- 



(3) 

The first statement is proven by observing that TfioodMEO does 
not depend on r ctr , therefore, as r ctl — > + , Pmeo converges 
to the optimal value of the centralized case. 

In order to prove the second statement, observe that 
diam(p[ 1 ](0),...,pN( )) < (n - l)r cmm and T FloodMEO G 
0(n). The result follows by substituting these bounds in (|3jl. 

■ 

Remark 5.4: The previous theorem confirms the intuitive 
idea that, if the communication is much faster than the mo- 
tion (r ctr small), then the optimal solution in the distributed 
case converges to the one of the centralized case. □ 



B. Distributed minimum time rendezvous in one dimension 



-6 



In one dimension (all the agents spread on a line), we can 
find a condition on r ctl ensuring that the move-toward-MBC 
algorithm is the solution of A4TTZ(Ea^ 1 B(Q, r ctr )). 

Theorem 5.5: For d = 1, let i max and i m j n the agents in 
the network S with the maximum and minimum positions. If 



r ct r < 



4 ' cmm 



and both max 



{p 1wJ (o)-p w (o)} > 

2r ctr and max jeA r(i mi]1 ){p l ' 7j ( ) ~ P 1 *'™ 1 (°)} ^ 2r ctr, then the 
control and communication law CCmeb solves the task 7^ n d zvs 
in minimum time which is exactly the time of the centralized 



solution, that is, T* = 



diam(p^ ,...,p^ n 'J 



□ 



Proof: Consider the input sequence of the centralized 
solution for the agent i m j n (and equivalently for i max ). It 
is u^(t) = r ctr for all t < T* - 1 and U ^{T* - 
1) = MBC(pW(0), . . . ,p [n] (0)) - P M(T* - 1). Since the 
rendezvous time is bounded by the time that i max and i m \ n 
take to reach MBC(pW(0), . . . ,p [n l(0)), we need to prove 
that, as long as the consensus on the minimal enclosing ball 
is not reached, then y} lmi "\ (t) = — ul 1 '™] (t) — r ctr . Due to 
the symmetry of the problem we will give the proof only 
for imin. It can be easily shown that for all t > 1 such 
that Pmal{t) 7^ p[ lm "l(0) (consensus is not reached), the 
following holds: 

Pmax(t + 1) > Pm'ax(t ~~ -0 + r cmm- 

It follows: 

MBC^(t + l) = ^(t + lJ+pW^), 

then 

MBCl*-](* + 1) > ^(pW(< - 1) + r cmm + P W(0)) 
= MBC [t ™" 1 (f-l) + ir cmm 
>MBC |l - ] (f)-r ctr +Jr cram . 



This leads to 



' ctr i 2 cmm 



> TW, and r ctr < -r c 



The other two assumptions ensure the condition for t = 0. 



C. Simulations 

In order to illustrate the performance of our rendezvous 
algorithms, we implemented the move-toward-MBC algo- 
rithm, based on the FloodMEB consensus algorithm. We 
implemented it in the plane, d = 2, over the disk graph. The 
simulation run is illustrated in Figure^ The 32 agents have a 
bound on the control inputs r ctr = 0.1, and a communication 
radius r cmm = 3. The initial positions of the agents were 
randomly generated over the rectangle [—6,6] x [—3,3]. 

The FloodMEB law converges in five steps, while the 
rendezvous is achieved at T = 58. As it clearly appears 
in the figure, once the consensus on the minimal enclosing 
ball is reached, all the agents move toward the center. 
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Fig. 1. Evolution of the network (in filled blue) according to CCmeb with 
evolution of FloodMEB (green circles connected by dashed red line) 



VI. Conclusions 

We have presented some simple algorithms on how to 
compute optimal enclosing shapes for pointsets via dis- 
tributed computation. These algorithms are then used to pro- 
vide efficient solutions to distributed rendezvous problems 
for synchronous robotic networks. For future work we en- 
vision characterizing the time complexity of the FloodMEB 
algorithm and, in turn, of the move-toward-MBC control and 
communication law. 
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